Helically symmetric iV-particle solutions in scalar gravity 
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Within a scalar model theory of gravity, where the interaction between particles is given by the 
half-retarded + half-advanced solution of the scalar wave equation, we consider an A^-body problem: 
we investigate configurations of N particles which form an equilateral iV-angle and are in helical 
motion about their common center. We prove that there exists a unique equilibrium configuration 
and compute the equilibrium radius explicitly in a post-Newtonian expansion. 
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Self-gravitating systems with helical symmetry have 
recently attracted considerable interest in the context of 
the numerical evolution of coalescing neutron star and 
black hole binaries 1, 2]. Since numerical codes simulat- 
ing relativistic collapse cannot evolve for large times, ini- 
tial data should be close to the final plunge. A reasonable 
approximation to such initial data might be to consider 
data given by spacetimes with a helical Killing vector, the 
reason being that gravitational radiation tends to reduce 
the eccentricity of orbits. However, not only for numer- 
ical purposes, but also from a more systematic point of 
view, spacetimes with a helical Killing vector form a class 
of time independent solutions of the Einstein equations 
which are interesting in themselves and about which lit- 
tle — including their existence — is known. The simplest 
class of examples are N point particles of equal mass m 
in Newtonian theory forming an equilateral A^-angle and 
uniformly rotating about their common center. We con- 
struct the analog of these solutions in a special relativistic 
scalar theory of gravity, with particles interacting via the 
half-retarded + half-advanced (" symmetric" ) solution of 
the wave equation. The case N = 2 (and allowing for 
different masses) has been considered by 0, 0| i the elec- 
tromagnetic N = 2 case has been treated in the seminal 
paper |5j. In this paper we perform a careful study of 
the symmetric interaction of particles in helical motion 
which is absent in the literature even in the antipodal 
(N = 2) case. This serves the purpose of proving the fol- 
lowing result: given m and the angular velocity Q of the 
helical motion, there exists, like in Newtonian theory, ex- 
actly one radius f c , for which the symmetric interaction 
is balanced by the centrifugal force. 

Although scalar theories of gravity are known to dis- 
agree with experiment, they provide simple test models 
for relativistic gravity, mainly because they have just one 
degree of freedom, see [6[ . Scalar theories derive from an 
action 

S = \ I g al3 $, a <Z>.^x + 4ir f P F{<5>)d A x, 

1 JM J M 

where (M, g a p) is Minkowski space and p the energy den- 
sity of matter; we have set G = 1, c = 1. The resulting 
equation for $ is □$ = 4npF'(&). When F(0) is chosen 



to be one, the associated theory has the correct Newto- 
nian limit. The choice F($>) — exp_$ corresponds to the 
model theory recently proposed in [6] . To further simplify 
matters, we make the choice = 1 + <f>, which corre- 

sponds to a first-order expansion, see [7|, |8[ ; this leads to 
the linear wave equation for $, 



□$ = 47rp . 



(1) 



We consider a family of N structureless point particles 
of equal mass m; let x n (s n ) be the world line of the n th 
particle, where s n denotes proper time; then 

AT-l „ 

p(x) = m / <5 (4) (x - x n (s n )) ds n , 

n=0 ^ 

so that the particle equations of motion are 

| Sii =0 (2) 



ds n 



for n = 0, 1, . . . (N— 1), where the field <i>| Sn acting on the 
n th particle is the symmetric solution of |T]) generated by 
the remaining particles 

In helical symmetry, fields are invariant under the ac- 
tion of a helical Killing vector £, whose components are 
£t = 1; £0 = ^ = const, £ r = = £ z , and whose Lorentz 
norm is £ 2 = —1 + fl 2 r 2 ; here, (r, <f>,z) are cylindrical 
coordinates associated with x = (x,y,z). Helical mo- 
tion is motion tangent to the Killing orbits, i.e., circu- 
lar motion with constant angular velocity £1 in planes 
z = const. When we define fj, = <f> — fit, we find that a 
field ip on (M,g a p) is helically symmetric, if it is of the 
form ipip-i r i z )j where ip is periodic in fi with period 27r; 
helical motion is any motion with (fx,r, z)(s) = const. 

Helical solutions of the wave equation. We consider 
the wave equation |T]) for a helically symmetric source 
p(t,x) = ph([i,r,z), where ph is 27r-periodic in ju,; we 
assume ph = for r > f2 _1 so that the source is con- 
fined within the light cylinder r = where velocities 
are less than the speed of light. The retarded solution 
& le t(t,x) and the advanced solution $ a dv(£, %) of {J) are 
given by 

p(tT \x-x'\,x') 



ret / adv 



(3) 
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The solutions $ re t and <j> a dv share the symmetries of the 
source, i.e., <& rct = & Ict ((J,,r, z). To make this explicit, 
we first introduce cylindrical coordinates (t, r, cf>, z) and 
(r',(j)',z') associated with (t,x) and x' in ([3]); we then 
find that the integrand can be regarded as a 27r-periodic 
function of a = d> — 6' which further entails that 



2-n fT 

$ ret = - j da I r'dr' / dz 



, p h (p - a + il\x - x'\,r' , z') 



\x — X' 



where \x — x'\ 2 — r 2 + r' 2 — 2rr' cos a + (z — z') 2 . Con- 
sequently, $ re t = $ret(M; r ; z ) with 27r-periodicity in p. 
Note that the advanced solution < &adv(M! r i z ) arises from 
'I'ret (m, r, z) by making the replacement f2 — > (—0). 

We proceed by defining a variable p' via 

2r/coscr+(z-z') 2 ]^. (4) 



<7 + n [r 2 +r' 2 



For fixed p, r, r' < f2 , z, z' , the map a i— ► // is mono- 
tonically decreasing and thus a diffeomorphism. This 
follows from a straightforward computation, where we 
invoke de l'Hospital's rule (for r = r' , z = z', a = 2kn, 
k G Z). Performing a change of integration variables 
from er to //, where we use that a shift by 2n in er causes 
a shift by — 27r in p! , we eventually arrive at 



Z7T \l OO 

$ rct = J dp' J r'dr' J dz' 



p h ((i,',r',z') 



p — p' — o~ + Q 2 rr' sin a 



In this integral, a is to be regarded as a function of the 
other variables, implicitly given by ([4]). In fact, a = 
a(p — //, r, r', z — z'), where er(/i, r, r', z) satisfies 



(j, — a 



9. 



(r - r') 2 + 4rr' sin 2 — + z A 



= 0. 



(5) 



We call a(p,r,r' , z) the retarded angle associated with 
p (and the particular choice of r, r' , z). The following 
properties of a are immediate from the above discussion: 

^>0, cr(p + 2tt, r, r' , z) = cr(p, r, r' , z) + 2tt . (6) 
dp 

Finally, regarding the retarded solution $ rot (/i, r, z) as 
the convolution of an integration kernel with the source 
Ph(fJ>' ,r' , z') yields the so-called retarded kernel 



K Tet (iJ.,r,r' \z) = Q, 



1 



a + fl 2 rr' sin a 



(7) 



where a = o~(fi,r,r' , z). As a consequence of (|6l). 
K let (p, r, r', z) is 27r-periodic in /i. 

As noted above, the advanced kernel -Kadv(Mi r i r 'i z ) is 
given in analogy to ([7]), where — > (— 0) and a — ► <7 a d v ; 



M — °adv 



(r - r') 2 + 4rr' sin 



. 2 ^ a dv 



From ([5]) we conclude that [— er(— /i, r, r', z)] satisfies this 
equation, hence <Jad v {p,,r,r' ,z) = —<r(—fi,r,r',z). By 
periodicity of K re t we thus infer the important relation 



K^ip^r^r'^z) = K rct (2ir - p,r,r',z) 



(8) 



Symmetric solution for a point source. The simplest 
source that is compatible with helical symmetry is a point 
mass m in circular motion. Let (fi,f < i7 _1 ,z = 0) be 
the position of the point particle; then the density is 

p h (/x, r,z) = m (1 - Q 2 f 2 ) h 5(p - p,)^jls(z) , 
and the associated retarded potential reads 

§ rct {fi,r,z) =m(l- Q 2 f 2 ) 2 K Ie t(n- p,r,f,z) . 

The radial component of the force at (ji, r, z) ^ (ji, f , z) 
is given by <9 r $ rc t; when r = f, z = it simplifies to 



(i-r» 2 f 2 ) 5 ^jr rot (^-/i > f,f,o), 



where we have used the symmetry of the retarded kernel 
in r and f, i.e., K Tet (fj,,r,f, z) = K Tet (fJ,,f,r, z). Con- 
sequently, the fields $ ro t and [9 r $ rc t]| r =r at positions 
jj, 7^ ft, r = f, z = are completely described by 
K Te t{fx, r, f, 0) and its derivatives. 

To obtain the kernel K rct (n, f, f, 0) we first investigate 
the retarded angle a(p,f,f,0). It ensues from ([5]) that 
/it = corresponds to a — and thus ^ = 2n to <r = 27r 
by ©. Hence, for /i G [0, 27r), sin ^ is non-negative and 
the defining equation for <r(/i, f , f , 0) thus becomes 



fx — a 



2Qr sin ■ 



. 



(9) 



We find tr = /i in the limit f — > 0; for f > 0, however, 
a > fi, unless cr = = /i or it = 2ir = ji; the difference 
between the angle fi and the retarded angle a is largest 
when a — n. Keeping \l e (0, 27r) fixed (so that a is in 
the same interval) we obtain 



dfCr 



2Q sin • 



1 - ftfcosf 



>0; 



(10) 



hence, in addition to being increasing in /i, a(fj,, f, f, 0) is 
increasing, with range (/x, o~i(p) < 27r), also when viewed 
as a function of f. 

Finally, the results © and (jTO]) lead to 
1 1 



2fsinf 1 



Qr cos ; 



11) 



8r Krr.t — 



(1 - Of) 2 +4fifsin 2 f 
2f 2 sin % 



1 



1 — Qr cos : 



3 



p! = 2ir/N 




FIG. 1: N point particles of equal mass m in uniform circular 
motion — the mutual interaction is given by the symmetric 
potential (half-retarded + half-advanced potential) . We prove 
that there exists a unique radius r e such that the configuration 
is in equilibrium. 



The kernel K rct is manifestly negative and monotonically 
increasing in r, since dfK let is positive; for f — ► both 
expressions diverge, while the limit is finite for f — > Vt . 

Since K rst (TT, r, f, 0) = #adv(7r, f, f, 0) by ©, the re- 
tarded and advanced potential at the antipodal point 
(fj, = p, + 7T, r = F, z = 0) are equal; the same is true 
for the radial components of the forces. The tangential 
components, however, are not equal, but opposite, since 
9/*ifretU=7r = -0/*ifadvU=7r by ©. Therefore, if we con- 
sider the symmetric solution $(/i,r, z) of i.e., 



$(^,r,z) 



$ ret (At, r, z) + $adv(A*, r, z) 



(12) 



then the /i-derivatives of the two terms cancel, so that 
= at the antipodal point of a point mass. This 
fact is a necessary prerequisite for a system of two (or 
more) particles in circular motion to be in equilibrium. 
Henceforth we only consider symmetric potentials (fT2|) . 



Equilibrium configuration for N point masses. We 
now consider the helical configuration depicted in Fig. [1] 
let n — 0, 1, . . . , (N — 1) be TV point masses of equal mass 
m, equidistantly distributed along a circle of radius f at 
z = and uniformly rotating about their common center 
- the n th particle's position is thus given by (/2„,F, 0) 
with /2„ = 2%n/N. Let $„(/i, r, z) denote the symmetric 
potential generated by the n th particle. At the position 
(fx, r, z) = (0, f, 0) of the first point mass the total poten- 
tial $ is then given as Yl n >x ®n{®, r, 0), hence 



$(0,r, 0) 



(1 - Q 2 f 2 ) 



JV-l 

£ 

71 = 1 



(i^rct +iC a d v )(/i n ,r, r, 0). 



Making use of ((H) results in 

AT-l 

$(0,r,0) =m(l-n 2 f 2 y K rct (fLn,f,f } 0) 

n=l 
N—l 

[d r ^]\ r = ^(l-il 2 r 2 yJ2 d r K ^(fin,f,f,0) (13) 

n=l 

for the potential and the radial component of the 
force at (0, f, 0). The tangential component of the 
force, i.e., d^(p, f, 0)| M=0 , vanishes, since 9 M A" a dvU = 



-^i^retlavr-zj and thus c^ifadvUn 



likewise, d z $(0, f, z)\ z= o — 0, which is a simple conse- 
quence of the mirror symmetry in z. The equation of 
motion ([2]) for the first particle thus reduces to 



[a r $(o,r,o)]| r=F 



fl + $ 



(0,f,0)) T 



2 r 



ft 2 r 2 



0, (14) 



where we have used the independence of (proper) time 
of $ at the particle's position and tuxq — (<9 Q £ 2 )/(2£ 2 ), 
where £ 2 = — 1 + £l 2 r 2 , for the centrifugal term. The 
equations for the remaining (N—l) particles are identical, 
since the symmetry of the configuration entails that none 
of the particles is distinguished. Hereby, the system of 
3iV equations © reduces to one single equation (fT4"]) . 

We conclude that the configuration of Fig. [T] is in equi- 
librium, if condition Q14p holds, i.e., if the radial force 
acting on each particle is balanced by the centrifugal 
force. As follows from (|11[) . the potential <&(0,r, 0) is 
negative and monotonically increasing for all f2r < 1; it 
diverges for f — ► and converges to zero when Qf — > 1. 
Consequently, there exists a unique radius fo such that 
(1 + <&(0,f, 0)) is negative for all f < fo, and positive 
for r > f . The term \d r $(0, r, 0)1 \ r= f is positive for all 
VLf < 1; it diverges as f — > and it goes to zero when 
ilf — ► 1. Combining the results it follows that the func- 
tion on the l.h.s. of (fTl)) is positive for all f < f and 
goes to — oo as Clf — > 1. We thus conclude that this func- 
tion assumes the value zero at least once in the interval 
f <E (fo,^ -1 ), so that there exists at least one radius f e 
for which condition (fT4)) is satisfied and the configuration 
is in equilibrium. In the following we prove that radius r e 
is unique by showing that the l.h.s. of (fl4|) is decreasing 
for f e (fo, fi™ 1 ). 

The proof would be trivial if the radial force 
[d r <&(0, r, 0)1 | r= f were decreasing in f (since the sec- 
ond term on the l.h.s. of fl4|) is manifestly decreas- 
ing for f > fo). However, whether monotonicity of 
\d r $>(0, r, 0)1 | r=f actually holds, is unclear in general. 
Namely, it can be shown numerically that VKMj ^) — 
df[(l — il 2 f 2 ) 1 / 2 9 F iir ro t(/i, r, f, 0)] does not have a sign: 
there exists a connected domain D in the set (0, 27t) x 
(0, SI -1 ) such that ip is positive when (/i,f) e D and 
negative when (ji, f) ^ D — this is in stark con- 
trast to the Newtonian case, where ip is negative for 



4 



all The main properties of D are the following: 

mine f w 3/4 f2 _1 , hence negativity holds for small f, 
where velocities are small compared to c so that Newto- 
nian gravity is a good approximation to scalar gravity; 
maxfl fi w 1/4, hence the radial force is decreasing at 
least when the number of particles is sufficiently small, 
i.e., when /2„ > 1/4 Vn ; typical values of ip on P> are 
by several orders of magnitude larger than typical values 
of on [(0, 2n) x (0, fi -1 )] \P> — this complicates mat- 
ters when one seeks to prove that the sum over all \x n 
is negative. Despite this last remark, numerical evidence 
suggests that [9 r $(0, r, 0)1 \ r -f is in fact decreasing ir- 
respective of the number of particles; a rigorous proof, 
however, seems difficult to obtain. 

In our proof we therefore proceed along different lines. 
The derivative of the function on the l.h.s. of (fT4l) reads 



7V-1 




df 



fiV y drK Tet (ji n ,f,f,o) 



) 2 K Tet (fl n ,f, f,0) 



n 2 f 



i - n 2 r 2 



(l + $(0,f,0)) d f f - 



o 2 



-n 2 f 2 



(15) 



Since the last line is clearly negative when f > Tq, in 
order to show that the whole function is negative, it suf- 
fices to prove that each of the terms in braces is negative 
individually. To this end let a = <r(fL n ,f, f, 0) for some 
n: then each individual term in braces has the form 



4f 3 (1 -ft 2 f 2 cos§) sin | 



P(flr, cos ■ 



(16) 



where P(i>, cos is a complicated polynomial of degree 
eight in v — Qf < 1 and of degree four in cos^-. In 
a second step we replace cos ^ by a variable 6 defined 
through cos ^ = u _1 [l — (1 — v 2 )Sj. Since — 1 < cos § < 1, 
the permitted range of S is 



1 

2<r 



■V 1 — V 

Using 5 leads to a simple representation of P(v, cos ^) 
3 - 85 + <S 2 (5 - 4« 2 ) + S 3 {2 + 4« 2 ) + 2S 4 v 4 



(17) 



P 



(1-w 2 )- 



The roots of the polynomial P are explicitly given by 



lie outside of the admissible domain (fT7|) of the variables 
(v, 5). Since in addition P < for v — » and <5 = 1, it 
follows that P is negative everywhere on the admissible 
(u, <5)-domain, or, equivalcntly, 



V(f,cr) € (O,^ 1 ) x [0,2tt) . (19) 



P(fif, cos-) < 



With P < the expression fj 16[) is negative, which com- 
pletes the proof of the claim. 

Post- Newtonian expansion. For a given number of 
particles, the equilibrium radius f e of the configuration 
in Fig. [1] is a function of the angular velocity and the 
mass. This functional dependence cannot be made ex- 
plicit, since this would involve, among other things, an 
explicit knowledge of the retarded angle © ■ (For a two- 
particle system, f e can be given as a (non-explicit) func- 
tion of the orbital velocity f2f e , see [3j, which, of course, 
does not lead to an explicit solution for f e .) 

It is feasible, however, to analyze the equilibrium con- 
dition (|14[) by means of a post-Newtonian approximation 
scheme. With the support of a computer algebra pro- 
gram necessary manipulations can be done in a straight- 
forward way and we eventually obtain a post-Newtonian 
expansion of f c ; here, we merely state some results. 

Let u) be the angular velocity as measured in standard 
units, i.e., [u] = s _1 ; clearly, ui = Vic, where c is the 
speed of light; furthermore, let G be the gravitational 
constant and M = Nm the total mass of the iV-particle 
configuration of Fig. [TJ We define a quanitity R (with 
unit length) and a dimensionless quantity x according to 

R = (GMlu- 2 ) * , x= (GMluc' 3 ) * . 
In terms of R, x the post-Newtonian expansion of f e is: 



N 


r e 


2 


\R [14 


- a; 2 /12 - 7a; 4 /72 + 0(x 6 )] 


3 




- 7x 2 /72 - 0.065a; 4 + 0(x 6 )] 


10 5 


1.228P [l - 


f 0.273a; 2 - 0.465a; 4 + Q{x 6 )] 



For highly relativistic configurations numerical investi- 
gations indicate that f c cx fl^ 1 as il — > oo, so there does 
not exist an innermost circular orbit. 



6 (1 — 8 ± vA), where the discriminant A reads 



A 



-2 5 



[V2-1]) (<5+[V2 + l]) . (18) 



Evidently, A is non- negative if and only if S < — 1 — \/2 
or 8 G [V2 — 1, |]. As a consequence, the roots of P 
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